Glassy effects in the swelling/collapse dynamics 
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Abstract 

We investigate, using numerical simulations and analytical arguments, a simple one dimen- 
sional model for the swelling or the collapse of a closed polymer chain of size N, representing 
the dynamical evolution of a polymer in a O-solvent that is rapidly changed into a good solvent 
(swelling) or a bad solvent (collapse). In the case of swelling, the density profile for interme- 
diate times is parabolic and expands in space as i 1 / 3 , as predicted by a Flory-like continuum 
theory. The dynamics slows down after a time oc N 2 when the chain becomes stretched, and 
the polymer gets stuck in metastable 'zig-zag' configurations, from which it escapes through 
thermal activation. The size of the polymer in the final stages is found to grow as Vint. In the 
case of collapse, the chain very quickly (after a time of order unity) breaks up into clusters of 
monomers ('pearls'). The evolution of the chain then proceeds through a slow growth of the 
size of these metastable clusters, again evolving as the logarithm of time. We enumerate the 
total number of metastable states as a function of the extension of the chain, and deduce from 
this computation that the radius of the chain should decrease as l/ln(lnt). We compute the 
total number of metastable states with a given value of the energy, and find that the complexity 
is non zero for arbitrary low energies. We also obtain the distribution of cluster sizes, that we 
compare to simple 'cut-in-two' coalescence models. Finally, we determine the aging properties 
of the dynamical structure. The subaging behaviour that we find is attributed to the tail of the 
distribution at small cluster sizes, corresponding to anomalously 'fast' clusters (as compared to 
the average). We argue that this mechanism for subaging might hold in other slowly coarsening 
systems. 
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1 Introduction 



The problem of the dynamics of a single polymer chain has been addressed in many ways over the 
past, the Rouse model being the simplest theory 0. In the case where the polymer is not a Gaussian 
chain, but evolves either towards a swollen or collapsed state, the dynamics can also be studied using 
more or less elaborate scenarios, either in a phenomenological way (like the "sausage" or "necklace" 
pictures in the case of the collapse to a compact state @,@|), via numerics or in an analytical 
way |J, Although a very common experimental realization of such a dynamical process is the 
folding or unfolding of proteins as they reach (respectively) their native or denatured state, it has 
been so far very difficult to test the various theoretical speculations [§]. 

The dynamics of polymer melts, on the other hand, is much easier to study experimentally M 
and can be compared to some computer simulations ||10|| . This case exhibits many interesting glassy 
features, but looks quite different from what happens to a single polymer chain. However, one can 
learn a lot from these studies, since for a long enough polymer chain that collapses onto itself into 
a compact conformation, the local environment of a monomer inside a globule is very similar to the 
one in a melt. It has recently been speculated that a single homopolymeric |11| or heteropolymeric 
chain |12[] does exhibit some glassy dynamics during the collapse or folding, including aging, at least 
in an intermediate regime of time before complete folding. 

While it is hard to find a satisfactory picture of what happens microscopically as a real polymer 
chain collapses in three dimensions (either because of the length of computational time for the 
numerics, or because of the complexity of the analytical approaches involved), one can gain some 
intuition by looking at the simpler (but less physical) case of a one-dimensional polymer chain. 

In this paper, we investigate theoretically and via Monte-Carlo simulations the dynamical be- 
haviour of a single one- dimensional polymer, composed of identical monomers that interact via local 
two-body interactions. Our aim is to mimic what happens to a polymeric chain in a good or a bad 
solvent when it swells towards an expanded coil, or collapses onto a compact globule. We will present 
some numerical results for both cases, together with some approximate analytical calculations which 
rationalize these numerical findings. For example, we show that in several situations, the chain con- 
straint (which is difficult to take into account in a rigorous way) can be neglected, and the polymeric 
chain behaves (at least in some intermediate time region) as an assembly of independent particles 
ruled by the same type of dynamics. 

Finally we will discuss how our results could be compared to real systems easier to observe, like 
polymer melts [JTTJ . The model we present here is in some way a toy model for non-disordered systems 
where the aggregation interactions are most important, with or without elastic chain constraint (a 
realization of that being colloidal gels where glassy cluster patterns appear at high concentrations 

PI)- 



2 Model and simulations 

We consider a periodic string of N particles (monomers), each of which constrained to live on the 
sites of a one-dimensional lattice. The elasticity of the chain is enforced by the fact that successive 
monomers along the chain cannot be at a distance larger than m lattice spacings in space. In our 
numerical simulations, we have chosen m = 2. If one allows m = oo then the system is no more an 
elastic chain but an assembly of N independent particles. The polymer is constrained on a segment 



of L sites, defining a density of monomers p = N/L. This length L is actually the measure of the 
extension of the chain; as we study the dynamics of the chain, this length will evolve with time, and 
will be referred to as L(t). We assume that the interaction is short range, i.e., two monomers can 
only interact when they stand on the same site. (The influence of longer ranged interactions on the 
physical picture obtained below will be discussed in the conclusion). We will call M the number of 
sites occupied by more than one monomer, and is the number of monomers sitting on site %. The 
total energy E of the chain is therefore: 

E = 1^(^-1). (1) 
z i=i 

The interaction parameter v is positive for an excluded volume interaction, which leads to swelling. 
The ground state and stationary solution of the dynamics is then the fully extended chain, with an 
end-to-end distance R ~ N. For v < 0, the interaction is attractive, and we study the collapse of the 
chain. In this case, the corresponding lowest-lying energy state is a compact "cluster" of N particles 
sitting on a unique site. As will be seen in the following, this final state is actually never reached in 
observable times, as the dynamics dramatically slows down with time. 

We shall choose as our initial condition a simple random walk for the chain, corresponding to an 
equilibrium configuration in the non-interacting case v — 0. The typical initial size of the chain is 
therefore ~ \fN . 

During the Monte-Carlo simulation, we do not allow collective motion of monomers but only 
individual moves. At each step, a monomer is allowed to move to one of its neighbouring sites as 
long as the distance between its new position and the neighbouring monomers along the chain is < m. 
In addition a monomer move is accepted according to a Metropolis criterion: If 5E is the energy 
difference between the new configuration and the initial configuration, then this move is accepted 
with probability p = min[l, exp(— ^)] where T is the temperature. If one monomer on site i is 
moved to site j, the energy change SE is given by: 

v 

5E = - [(rij + l)rij + (rii - l)(rii - 2) - n^rij - 1) - 71,(71* - 1)] = v(rij - n* + 1). (2) 

In the following, we will often choose v — ±1, corresponding respectively to swelling or to collapse. 
The quantities of interest are primarily based on the density profile n(x,t) = rii(t), where x = ia (a 
is the lattice spacing). For example, one can study the participation ratios Y q (t), defined as: 

X 

Obviously, Y\ — N and E = v(Yi — Y\)j1. Yq can be taken as a measure of the size of the chain. 
Alternatively, we will consider the gyration radius R(t) defined as: 



N , 



Note that for a fully stretched chain, which is the ground state of the chain, the maximum value of 
R is, for m = 2, N/Vl2 ~ 0.28N. 



3 Repulsive case: swelling of the chain 



3.1 Numerical results 

We have investigated the behaviour of the gyration radius R(t) for different chain sizes and different 
temperatures T, including T = 0. Interestingly, the different results can be rescaled on top of each 
other by plotting r = R/N as a function of r = t/N 2 : see Fig. I. There is a well defined regime 
of time scales where one finds r oc r 1 ^, with z ~ 3, which holds when r < 1. When r becomes of 
order 1, the chain becomes substantially stretched, and the dynamics stops (at zero temperature) 
or becomes activated (for non zero temperatures). The scaling r oc r 1 / 3 does not hold at very short 
times either, but sets in after a time t c ~ VN, as can be seen from Fig. 2 where we have plotted 
[R(t) — R(t = 0)]/VN as a function of (t/y/N) for short times and different values of N. The short 
time behaviour is linear in time. This initial time regime is found to be independent of the chain 
constraint, up to a rescaling of time. This is expected from the simple theories developed below. The 
chain constraint only starts to play a role when r ~ 1, i.e. for t ~ N 2 . 

The dynamics of the whole density profile n(x, t) from the initial configuration is also interesting: 
after a short transient time t c (< 100 MC steps for N = fOOOO) the profile n(x,t) smoothes to a 
parabolic form that expands gradually in space (see Figure 3). This expansion is self-similar in time 
and can be well fitted at large times by: 

(3) 

where ±^y^ are the points where the density vanishes. We find, for iV = 50000, n (t) ~ 725 1~ - 3 
and L(t) ~ 56.8 t a35 , compatible with z = 3. This form suggests that n (t) oc 1/L(t), as expected 
for a scaling profile. The participation ratios indeed read in this case: 

n = j2^ nl{t)L{t) K t(1 ~ q)/3 - (4) 

In particular, Y 1 = N = 2n (t)L(t)/3, whereas the energy is expected to relax as t~ 1//3 , a result that 
can be checked numerically independently. For such a parabolic profile, the mean-squared radius is 
given by R = L(t)/2y/E. 

Again, it is important to note that, except for the very last stages of the swelling, the results 
are very similar whether the chain constraint is present or not; we will discuss this feature also in 
the following section. The final regime towards the completely expanded configuration of the chain 
is hard to study numerically as the chain gets trapped in some metastable states. The relaxation 
towards the completely swollen state is slower and slower, as can be seen on Figure f . 

In the next section, we present different analytical approaches to describe the swelling mechanism 
to account for the above numerical results. 

3.2 Analytical theories for the swelling dynamics 
3.2.1 A Flory theory 

A very crude analytical approach for the dynamics of the chain can be derived from a dissipative 
equation of motion for the size of the polymer R, combined with a Flory approximation for the free 
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Figure 1: In this graph we plot r = R/N versus r = t/N 2 for various chain lengths and various 
temperatures. These are all single runs (no averaging over initial conditions). One can see that there 
is a well characterized r 1//3 regime for intermediate times. At zero temperature, the chain gets stuck 
in a metastable state and r saturates to a value much smaller than its maximum value r max = 0.28 
(shown as a dashed line). At non zero temperature, the long time regime is logarithmic in time. 
We have also shown, for comparison, the r 1//14 behaviour predicted by the dynamical effective length 
method discussed below in section 3.2.2. 
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Figure 2: Initial time regime, where (R — Ro)/y/~N is plotted as a function of t/yN for different 
values of N. One sees a first regime linear in t crossing over for t ~ to the t 1 / 3 regime. The 
dot-dashed line shows a fit to the simple Flory theory described below, Eq. (Q). Note that the curves 
are single runs, which explains the differences for short times. 




Figure 3: Profiles n(x,t) of the chain as it swells from a random initial configuration, at times 
t — 1, 10, 10 2 , 10 5 , with N = 10000. The profile quickly becomes parabolic and expands in space 
as t 1 / 3 . 



energy J r (R). The Flory free energy is as usual written as the sum of an entropic contribution and 
an excluded volume interaction contribution. The resulting equation reads: 



dR dJ r (R) 

f^chain' 



dt 



dR 



N dR 



R 2 N 2 



(5) 



where [i c hain is the mobility of the chain. It is well known || that the mobility for a chain is 
^chain = , where /io is the mobility of a single monomer. The exact solution for this differential 
equation can actually be computed, and for an initial condition R(t = 0) = Ro, reads: 



R(t) 



n' 1 



Ri +2 + 



vdN 3 



2(d+2) M0 t 

(e n 2 — 1 



1 

d+2 



(6) 



For large times, R(t) converges exponentially fast towards the equilibrium Flory radius Rp 

(vdN 3 /2)^: 



R(t) ~ R F 



2(d+2) M0 i 



d + 2 



(7) 



For an ideal chain initial condition, one has Ro ~ \HSf . Therefore, in the regime where N d ^ 2 /v <^ 
[i$t A^ 2 , the initial condition term R(0) d+2 can be neglected, and one finds: 



(did + 2) V+ 2 
R(t) ~ — -vfi Nt 



Note that in this time regime, the elastic part in the Flory free energy (that comes from the chain 
constraint) does not play any role and one therefore expects the same scaling to be valid for an 
assembly of non-connected repulsive particles, in the intermediate time regime mentioned above. As 
discussed above, this is indeed what we find numerically, up to a rescaling of time. 

In one dimension, the Flory theory predicts a power law behavior R(t) oc A^3t3 - both scalings 
in rather good agreement with the numerical results. The very short time behaviour predicted by 
this approach reads: 

*«^(0)(l + |p) (9) 

In one dimension, the departure from the initial condition is therefore R(t) = R(0) + r^t 7 ^ where 
jf = 1 and is independent of the number of monomers N. This approach therefore predicts that 
for short times R(t)/y/N is a function of t/y/N, in reasonable agreement with our numerical data 
(see Fig. 2). 

A slightly more refined theory (but still in the spirit of Flory) actually allows one to understand 
the parabolic shape of the density profile. One can write a Master equation for the number of 
particles sitting on site i: 

f)<n ■ 

^EWWO-EWW*) (io) 



where the transition probabilities at finite temperature T = 4 are chosen to be: 



In the zero temperature limit, the transition probabilities are simply described by a step function: 



Wji = W^j = w if Tii(t) > Ujit) + 1 

= if rii(t) < nj(t) + 1 (12) 

In the simplest case where we allow moves to both nearest-neighbour sites, we require the con- 
dition that Wij = Wji — if \i — j\ > 2. We neglect the chain condition, which is expected to be a 
valid approximation in the short time limit. The Master equation then becomes, in one dimension 
(d=l): 

dfi ■ 

= w n i+ i(t)Q(ni(t) - n i+ i(t)) + w ^-i(t)O(^(t) - n^t)) 

- w m(t) [&(n l+1 (t) - m(t)) + e(ni_i(t) - m(t))] (13) 

Let us assume that the density rii(t) is slowly varying in space, so that one can expand the function 
in powers of An = rii+i — rii as: 0(An) ~ 9 — /3vAnQi. Then, 

~ [n i+1 (t) + 7ii_i(t) - 2n,(t)] + ^[nf +1 (t) + n^^t) - 2n?(t)], (14) 

where we have rescaled the time by u>o©o, an d set (7 = fivQi/Qo. Taking the continuous limit both 
in space, we end up with the following non-linear conservative diffusion equation: 



dn(x, t) 
dt 



V 2 n + 2gV(nVn) . (15) 



Let us note that the assumption (3v\An\ <C 1 is crucial in the derivation of this continuous 
equation. This may be acceptable in the case of repulsive interactions (v > 0), where the profiles 
n(x, t) turn out to be smooth, so that the number of monomers on neighbouring sites is almost the 
same. Surprisingly, however, this non linear equation describes also very well the zero temperature 
results in the long time limit. This is not true in the case of attractive interactions (v < 0): the 
above equation then becomes unstable. This correspond to a 'pearling instability' in the collapse 
regime which we will discuss further in Section 4. 

The non-linear diffusion equation ( |I~5|) admits a self-similar solution for the profile n(x,t), of the 
form: 

^=wAm) (16) 

where R s (t) = t a N s is the scaling form for the end-to-end distance of the chain in the intermediate 
time regime that we study, and / a scaling function independent of both t and N. Introducing the 
rescaled variable u = x/R s (t), we find the following differential equation for u: 

- m^\fw + u f»] = + 2 3^mu) + /(«)/»] (i 7) 

We assume that in the regime of time studied here, the diffusion term can be neglected, which is 
true if t a < gN 1 ' 6 , and can be checked a posteriori when R s -C N. Then the equation reduces to: 

- aN'-H^ifiu) + uf{u)\ = 2gN 2 -' s [f'\u) + /(«)/»] (18) 



This fixes the values of both exponents: atp — g and 6f — 3, which gives back the results of the 
Flory equation of motion for d — 1. Furthermore, the resulting ordinary differential equation for 
f(u) has a parabolic solution. The final result, normalized to give the correct number of monomers, 
is Eq. (D, with n (t) = N 2 /' 3 t 1 / 3 (3/645-) l / 3 and L{t) = 2(9 g At) 1 / 3 , in agreement with the numerical 
results of Fig. 3 (obtained at zero temperature). 



3.2.2 A dynamical variational method 

An alternative, approximate method to study the dynamics of swelling and collapse of polymer 
chains was proposed in ||. It is based on a dynamical extension of the effective Kuhn length method 
of Edwards, which is known to reproduce the Flory value for the end-to-end distance exponent. 
The essence of the method, together with the self-consistent dynamical equations, are recalled in 
Appendix A. From these equations, the short time and long time dynamics of the chain was obtained 
in |J. In the short time limit, one finds that the chain radius grows, in one dimension, as: 

R(t) = R + T v t^ lv = -^ ocA°, (19) 

at variance both with the Flory prediction 7^ = 1, and with the numerical simulations. 

The intermediate self-similar expansion regime was not worked out in Looking for a self 
similar solution of the form R{t) oc t a N s , we have found (after some work) the following values for 
a and 5: 

a v = - 5 V = -, (20) 

Quite surprisingly, again, these values differ from the simple Flory prediction and are very far from 
the numerical result. From the above value of the exponents, one predicts a stretching time of the 
order of A 3 (longer than what is observed numerically, where the stretched regime is obtained when 
t = t/N 2 becomes of order 1). 

It is intriguing to see that a somewhat more refined theory (but which basically contains also 
the ingredients of a Flory theory) gives worse results than the simple-minded approach. It seems 
actually that the Gaussian variational method misses the fact that the dynamics is governed by 
repulsive interactions but not sensitive (at short times) to the chain constraint. Hence, the 'effective' 
chain approximation is very bad in this region. The approximation might be better at long times, 
to describe the approach to the swollen state when the chain constraint becomes relevant. We have 
shown in Fig. 1 the predicted t 1 ' 14 for comparison; however, as discussed in the next paragraph, a 
logarithmic growth is probably more adapted to describe this regime. 



3.3 Late stages of the dynamics 

It can be observed from the Monte Carlo simulations (see Figure 1) that the dynamical evolution 
becomes frozen at T = 0, while it is strongly slowed down at finite temperatures. In Figure 4, we 
show snapshots of the configuration of a short chain, with A = 100, as it evolves with time. 

It is clear that the chain is trapped for long period of times in metastable configurations, consisting 
of fully stretched segments of the chain going alternatively to the left or to the right. If we call 2M 
the total number of segments and I their average size, the typical extension of the chain in such a 
configuration is R ~ £a/M. Since A = 2M£, one has: R ~ N/\fM, or M ~ r~ 2 . The dynamics of 
the chain in the long time region therefore consists of a slow 'coarsening', where 'kinks' progressively 
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Figure 4: Snapshots of the configuration of a short chain, with N = 100, as it evolves with time. The 
x axis is the internal coordinate of the monomers, while the y axis is the corresponding position in 
space. One can see that the long time dynamics consists of the progressive disappearance of 'kinks'. 
Inset: The complexity S as a function of R. 



disappear, leading to a decrease of M. The time needed for a hairpin of size £ to disentangle can be 
estimated from the energy barrier to be crossed, which is proportional to £ itself. Therefore, one has 
lnt(£) ~ f3v£ ~ /3vNr 2 . From this argument, we obtain a logarithmic growth law of the size of the 
chain in the long time regime: 



R 



IN_ 
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hit, 



(21) 



in qualitative agreement with the numerical simulations: We show on Fig. 5 a logarithmic fit of the 
radius of gyration as the last stages of the stretching occur. The above estimate assumes that all 
segments lengths are close to the average value. This is confirmed by an analytical calculation of the 
number of metastable states of the stretching chain, and is explained in detail in Appendix B. 
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Figure 5: Logarithmic growth of the chain at the end of the swelling dynamics, for a small chain 
N = 100 



It is indeed interesting to count the number M{M) of metastable states of the chain for a fixed 
value of M, i.e., for a given size R of the chain. The problem is actually related to counting the 
metastable states of the one dimensional Ising model with Kawasaki dynamics; the corresponding 
calculations are given in Appendix B. The idea is to make an analogy between a segment of the chain 
going to the right (or to the left) and a domain of plus (or minus) spins. An additional constraint 
is needed to make the configuration metastable. The chain can only evolve by moving inwards 
the monomer at the extremity of the hairpin (see Fig. 4). This is possible if the total number of 



monomers sitting on the target site is less than the number of monomers at the starting site. If we 
assume that the local density of monomers is random, the probability for this to happen is equal to 
1/2. Therefore, a zig-zag configuration with M segments is stable with probability 2~ M . The result 
is that M{M) ~ exp(NS(M)), with a complexity S given by: 

S = -51n2 + (5-l)\n(l-5)-51n5, (22) 

with 5 = 2M/N ~ N/R 2 . This function is plotted in the inset of Figure 4. At the beginning of the 
swelling, the complexity increases very rapidly with R ; after a maximum reached for 5* = 1/3, S 
then decreases at large R. 

If one adds all frozen configurations for all possible values of M (or equivalently, R), one finds 
a total number of frozen states which grows as C N / \/N where the complexity per monomer is 
C = 3/2. This must be compared to the total number of configurations, which is equal, for m = 1, 
to the number of closed random walk, i.e. 2 N / y/N. 



4 Attractive case: collapse of the chain 
4.1 Numerical Results 

The case of attractive interactions can be studied numerically in the same way as before, except that 
now the interaction parameter v is negative. In this case, obviously, the ground state of the chain is 
such that all monomers are located on the same site. As the chain collapses, the profile n(x, t) very 
quickly becomes discontinuous as monomers aggregate in well separated "clusters": see Figure 6. 

These clusters are very slow to coalesce and the dynamics is slower and slower with time (and 
as the temperature T decreases). This slowing down is due to the fact that it is energetically very 
unfavorable to 'break' one of these clusters in order to transfer its monomers to larger clusters. The 
energy needed to remove one monomer from a cluster containing n monomers is equal to n\v\ (for 
n large), corresponding to a time cx exp(/3\v\n) . Therefore, the system remains trapped longer and 
longer in configurations where there are M clusters containing n ~ N/M particles. This can be seen 
for example on the time evolution of the energy E(t): single runs for three different temperatures 
are shown in Figure 7. By averaging E(t) over several MC runs, we found a logarithmic decay with 
time, compatible with the estimate E ~ —Mn 2 ~ — iVln(t), using n ~ Tint. It is interesting to 
note that the dynamics of clusters is very close to that of the one dimensional 'Backgammon' model 
introduced [l£j and further studied in [21, 22, 2C]: although the barriers are 'entropic' in the latter 



model, the average time needed to empty a cluster also grows exponentially with its size n. 
4.2 Breakdown of the continuum description 

Let us come back to the continuous approach developed in the previous section in the case of repulsive 
interactions. Starting again from the master equation, and considering attractive interactions this 
time, we find the following equation for n(x,t): 



dn(x, t) 
dt 



V 2 n-2\g\V(nVn) . (23) 




Figure 6: Profiles n(x,t) of the chain as it collapses from a random initial configuration, for T = 
at times t — \, 10, 10 5 . The profiles obtained for t = 10 2 , 10 3 and 10 4 are indistinguishable to the 
from the last one. 
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Figure 7: Energy of the polymer as a function of time for a single Monte-Carlo run, at T = 5, 50, 100: 
energetic traps can be visualized. 



A naive extension of the self-similar solution found in the repulsive case leads to the following equation 
for the size of the polymer R s : 

^ = -2M— (24) 

dt m R 2 s (ty [ ] 

or Rg(t) = R 3 (0) — QN\g\t, suggesting that all monomers coalesce into a single cluster in a finite time. 
This clearly does not capture the essential physical feature of the dynamics at late time, namely the 
activated coarsening of well separated 'towers' of monomers. 

In fact, Eq. (|23|) leads to a short wavelength 'pearling' instability of initially smooth profiles. 
Starting from n(x, t) = n + Sn(x, t), with 5n <C n , one finds: 

^^l = (l-2\g\n )V 2 6n, (25) 

which is unstable as soon as 2|g|n > 1, i.e. for small enough temperatures or when the chain density 
is large. As observed numerically, a random walk initial condition is indeed immediately unstable 
for large N, since n ~ VW 3> 1. 

As soon as the chain has entered the pearled state with well separated clusters of monomers, 
any continuous description is bound to fail. A discussion in terms of metastable states, and jumps 
between these states, is needed. 

4.3 Calculation of the number of metastable states 

A metastable state of lifetime t is a configuration of clusters which remains frozen for a time t. 
This means that (i) two clusters cannot be nearest neighbours and (ii) the minimum number of 
monomers in each of the clusters at temperature T is at least equal to k*(t,T) = max(2, Tln(t)) 
- since a cluster must contain at least two monomers. How many such configurations are there ? 
This is a combinatorial problem including two steps: first count the number of configurations with N 
monomers and M clusters each containing more than k* monomers, then the number of arrangements 
of M clusters on a one-dimensional box of L sites, so that all clusters are at least separated by an 
empty site. We will restrict to the case where the chain constraint is m = 1, i.e., two consecutive 
monomers cannot be at a distance larger than 1. In this case, the number of monomers within a 
'hole' is exactly equal to 2, because the chain is closed onto itself. Therefore, the total number of 
monomers belonging to clusters is related to the length L occupied by the chain by: 

M 

Y,n { = N -2(L- M). (26) 
i=i 

The number of configurations A4* (M\N) with at least k* monomers in each cluster can be calculated 
using the following formula: 

N N / M \ 

M k *(M\N)= E ••• E s(j2{n i -2)-N + 2L\, (27) 

n\=k* riM=k* \i=l / 

which can be transformed into: 

/ N \^ / i(k*—2)\\M 

A4* (M\N) = J d\e-^ N - 2L ^ I E e iA(n " 2) I = J d\ e -* XN I & - _ — x 1 (28) 



We define in the following 5 = 41, p = N/L and z = iX, and perform a saddle point approximation 
on the above integrand in the large N limit. We find that at the saddle point, z* = In 
Then, we obtain: 



l-2-<5(fc*_2) 
l-|-<5(fc*-3) 



ln[A4*(M|AT)] 



-l + - + <y(jfe* - 2) 
p 



l-Z-6(k*-2)\ ( 5 

Ln 1 1 _ i _ g (fc . _ 3 ) ) ~ Jln ( i - a - g(fc» - 3) 1 ■ (29) 



Now we count the number of ways A/"l(M) to arrange M clusters on L sites with all the clusters 
at least separated by an empty site. This can be done in the following manner, mapping this problem 
into an antiferromagnetic spin representation. (A more direct combinatorial method can be found in 
p0| ). We define the formal antiferromagnetic hamiltonian: 

L 

H = Jj2wi+i, J >°> (3°) 

where Oi = 1 by definition if there is at least one monomer on site i, and <7j = if site i is empty, so 
that configurations with consecutive occupied sites have a positive energy and the configurations of 
zero energy are those with at least one empty site between each occupied site. Using the additional 
constraint that the total number of occupied sites is set to some value M, the number of configurations 
7Vl(M) is given by the following formula: 

Al(M) = lim J2 e~ w 5(M -$>*)• (31) 

J— >oo , — , — * 

K=o,i} i=1 
This can be rewritten using the integral representation of the 5-function: 

Af L (M) = J d^ M Z L (ip) (32) 

where Z L {ip) = Y^{ ai =o,i} e~ w e~ v 2i=i°\ By performing a saddle-point approximation on /1, one has 
to find /i* such that 

M+-^| M=M *lnZ r (i A t) = 0. (33) 

The quantity Zi{ip) can actually be computed exactly using a transfer- matrix method, and in the 
limit J — > 00, Z L (ip) = s(l + vl + 4e _v . Using the saddle-point equation, /1* is easily found to 
be equal to vl-f 4e~ 4At = 1-2|5 . Finally, 

A/" L (M) ~e LS(<5p) , (34) 

with 

*W ^ **T + ™ ^ (irw) + k (l^) ' (35) 

One can check in particular that one finds S(0) = and £(§) = 0, as it should. 



We now combine the two above results, and deduce the total number of fc*-stable configurations 
for a given density p: 

N 

Af p (k*) = e NS "^ = Y, N L (M)M k *(M\N). (36) 

M=l 

This is again computed from a saddle point (in M) and leads to: 

2 (l-Z-6(k*-2)\ ( 5 

Sp{n = (-l + - + ^--2))ln( 1 _;_^._3 ) )-^ln ^_,_^,_ 3 ) 

- ' (37) 

In this expression, 5* is the saddle point value of 5, defined as: 

(1 _ 2 _ _ 2 ))**-*(l - 2 _ 5 *(k* - 3)) 3 - fc * (i _ 2 <Tp) 2 1 



(3* 



At zero temperature, all metastable states are frozen, and the size of the chain L and therefore the 
density p, are constant in time. One can estimate the complexity by setting k* = 2 in the previous 
formulae. We found numerically the solutions of the above equations 5*(p) and S(p) and plotted 
the results in Figure 8. One should note that since each cluster contains at least 2 monomers, the 
maximum number M of clusters is y, so S < |. Since a metastable state is such that two clusters 
should be separated by at least one empty site, one also has L > 2M, or 5p = ^ < |. 

It is interesting to draw an analogy with some recent work on the number of metastable states in 
spin-glasses Q. As k* increases, one counts configurations which are more and more stable. Similarly, 
in spin-glasses, one can introduce fc*-spin flip stable configurations, i.e. configurations the energy 
of which cannot be lowered by any flip involving at most k* spins. The complexity of these states 
was determined for simple spin-glasses by Biroli and Monasson They show that in mean-field 
models, inherent states (corresponding to 1-spin flip stable states), pure states and k* stable states 
all coincide. But they also show that in finite dimensions (contrarily to infinite connectivity models) 
all these different concepts are distinct and actually help to visualize the degree of freezing of a glassy 
system, according to the lifetime of the metastable states relevant after a certain time scale. In the 
following section, we will show that the analysis of the relevant metastable states of a collapsing 
polymer allows one to predict how the radius of gyration of the chain decreases with time. 



4.4 Consequences for the evolution of the chain at non zero tempera- 
tures 

At non zero temperatures, the relevant metastable states after a time t are such that k* ~ Tint. 
Furthermore, it is reasonable to assume that the chain size is such that the number of metastable 
states is maximal, that is to say, the chain is in one of its most probable configurations after a time 
t. We thus maximize the entropy and obtain p(t) = p*(t) (and therefore L(t) = N/p*(t)) by setting: 

dS p {k*) 
dp 



n=rr 



(39) 




Figure 8: Entropy of metastable states S(p) and the corresponding optimal number of clusters S*(p), 
as a function of the density of monomers p, at T = 0. 



Solving for p* and 5* in the large u = Tint limit, we find: 



p*(w)~41n-u ~ — 



1 



■u w lnw 

and the corresponding value of the entropy, to leading order in u: 

, „lnw ln(lnw) 1 , , , 
5* (u) ~ 2 2^ '- + - 4 - In 8). 

u u u 



(40) 



(41) 



These results predict a slow decrease of the number of clusters and an extremely slow decrease 
of the size of the chain: 



M(t) 



Tint 



L(t) 



41n(Tlnt) 



S(t) 



21n(T lnt) 
Tint 



(42) 



This last result is quite interesting and is a non trivial consequence of the present calculation of 
the complexity of metastable states. It is obviously very hard to check numerically, but suggests 
an hyper-slow contraction of the chain. Numerically, we find that the chain length indeed hardly 
decreases with time. 
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Figure 9: Entropy of metastable states S*(u) as a function of u — Tm(t). Inset: Density of the chain 
as a function of ln(-u) (plain line). The dotted line is the asymptotic result p{u) ~ 41n(w). 



4.5 Calculation of the density of metastable states at fixed energy 



We now turn to the calculation of the number of metastable states corresponding to a certain degree 
of freezing k*, and to a fixed total energy of the system E. As was done before in the case of 
spin-glasses [14] or Josephson arrays [15|, we want to understand how the entropy of metastable 



states behaves as a function of energy at T = 0, and we go here further by looking at its dynamical 
evolution with respect to k* = Tint. The calculation proceeds much in the same way as in the 
previous section (we also restrict ourselves to the case where m = 1), starting from the expression of 
the number of fc*-stable configurations with M clusters and N monomers, at fixed energy E: 



N N / M \ / M i 

M P ,AE,M\N)= £ ... £ <y £(„ i - 2 )-JV + 2L <f 

m=k* nja=k* \i=l / \i=l 



2 m(rH-l)-(L-M)-Ej . 

(43) 

The technical details of the calculation are relegated to Appendix C. Setting e = J| < and S p (e, k*) 
the complexity per monomer, our final result is plotted on Figure 10. We see that all S p (e, k*) for 
different k* follow a common envelope for a given p, and asymptotically the curves approach the 
following law at very low energies: 

S P (e) -^(1 - 2 -f + -(1 " -) 2 ["1 + |ln(f (1 - -) 3/2 )] + o(i) (44) 

e p e p 2 4 p e z 

For large times, p —>■ oo and we are left with: 

<? ( e) -^t^ (45) 

e 

Even for very low energies, the entropy of metastable states never reaches zero: contrarily to the 
example of spin-glasses, there is no lower cut-off in the energy per monomer below which the number 
of metastable states is not exponentially large in N. 



4.6 Distribution of cluster sizes 

Let us define p(n, t) as the distribution of sizes of the clusters, at a given time t and temperature, 
n being the number of monomers in a given metastable cluster. We further define n as the aver- 
age number of monomers in a cluster. From the above sections, we know that n(t) grows slowly 
(logarithmically) with time. We have studied np(n,t) as a function of z — n/n at different times. 
These different curves collapse well onto a single one (note however that the range of variation of 
n cannot be very large). The scaling function p is very far from a simple exponential, as a naive 
maximum entropy argument would suggest. In particular, the tail of p can be fitted as a power-law 
p(z) ~ (z + £o) 1 ~ u -i with v ~ 1.9. For z small, on the other hand, p(z) rapidly vanishes. Although 
our statistics is not very good, one can fit hxp{z) as (zi/z) c , with c ~ 1. 

The present cluster dynamics is actually quite similar to the so-called 'cut-in-two' model intro- 
duced by Derrida et al. ]17] to describe the pattern of coalescing droplets in one dimension. In 
the latter model, a collection of intervals of different sizes (which represent the droplets) evolves 
according to the following simple rule: one picks the smallest interval, divides it in two equal parts 
and sticks the right part to the right neighbour, while the left part coalesces with the left neighbour. 
In this model, the average size of the intervals grows with time, and the rescaled distribution of 




Figure 10: Entropy of metastable states S p (e) as a function of the energy per monomer. The plain 
line is the asymptotic law 5(e) = — ln(— e)/e. For a given k*, the domain of definition of the entropy 
is bounded by e max (k*) = — + 1) + — . The successive arches, defined on [e max (k* + 1), e max (k*)}, 
correspond from right to left, to increasing k* (= 3,4, ..), and follow a common envelope. 



interval sizes tends to an asymptotic distribution shown in Fig. 11. The cluster dynamics studied 
here is similar for the following reason: because the time needed to empty one cluster is exponential 
in its size, the first cluster to disappear will typically be the smallest one available after time t. If 
one neglects the chain constraint, the smallest cluster empties itself by sending an equal number 
n/2 (on average) of monomers towards its two immediate neighbours. The analogy is not exact, 
though, since this number is not exactly equal to n/2 at each iteration; furthermore, some subtle 
correlations are induced by the chain constraint. One can nevertheless take the simple 'cut-in-two' 
model, or a randomized version where a random fraction / is pasted to the left and 1 — / to the 
right, as benchmarks to which we can compare our results. This is performed in Fig. 10. The 
agreement is only fair, and is better for the randomized model than for the strict cut-in-two version, 
for which p(z) is known to decay exponentially for large z [17[]. Note however that the density of 



small clusters strictly vanishes below a certain value z c in the cut-in-two model, whereas it is finite 
in the present dynamics. This is due to the fact that although small clusters are typically the first 
ones to disappear, fluctuations can persist and keep some small clusters alive, whereas by definition 
they systematically disappear with a cut-in-two rule. The excess density of small clusters has the 
interesting consequence, discussed below, to give rise to a sub-aging behaviour. 
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Figure 11: Comparison of the distribution of the cluster sizes in our model with the ones correspond- 
ing to a strict cut-in-two model and a randomized cut-in-two model. Note that the density of small 
clusters strictly vanishes in the latter two models. 



4.7 Aging at long times 



Another manifestation of the long lived metastable states is the aging phenomenon, which we have 
investigated numerically by computing the following age dependent structure factor |[23|| : 

1 N 

C q (t w , t + t w ) = - £< e *(**(*+*«>-**(*«))) (46) 
^ i=i 

This quantity is expected to depend both on t w and t (violation of time translation invariance) in 
the out-of-equilibrium regime where T login, <C N. As can be seen from Figure 127, obtained for 
T = 5, N = 10000 and averaged over 10 initial configurations, the relaxation is slower and slower 



as t w increases. As often observed experimentally |24| , [2q , [26|1 , we have found that the data can be 
rescaled in the variable s = t/tfa, with \i < 1 (subaging). Within the time regime probed by the MC, 
\i ~ 0.75 for q = 5 and q = 10. This means that the characteristic relaxation time grows more slowly 
that the waiting time itself. In the present situation, this can be interpreted following the above 
remark about the small z behaviour of p(z): the presence of clusters much smaller than the average 
size n means that some clusters will evolve on a time scale much shorter than t w = exp(n/T). More 
precisely, one can approximate the short time behaviour of C q (t + t w ,t w ) by: 



1 - C q (t + t w , t w ) ~t J dz zp(z) exp(-zn/T). (47) 

The above equation is obtained by assuming that each cluster of size n contributes to the correlation 
function as nexp(— i/r(n)) where r(n) = exp(n/T). If hip(z) behaves as (z /z) c for small z, a 
saddle-point approximation of the above integral leads to: 

l-C^t + t^Uoc-, /x=(l + -) l^f) • (48) 

This scenario therefore leads to a subaging behaviour for large t w , albeit with a (slowly) time depen- 
dent exponent /i. The absence of 'small' clusters would correspond to the limit c — > oo, such that 
p(z < 1) = 0. In this case, one finds \i = 1 as expected. The rescaled function C q (s) can be fitted by 
C/(s + s) w , with w ~ 0.18 for q = 5 and w ~ 0.10 for q = 10. 

The above mechanism for subaging might be much more general, and hold for other slowly 
coarsening systems, such as the Random Field Ising model. Suppose that the distribution of domain 
sizes scales with an age dependent average length R(t w ), such that t w ~ exp{R^/T) 0. Then the 
presence of domains smaller than R(t w ) will correspond to relaxation times much shorter than t w , 
and therefore to the possibility of subaging effects. 



5 Summary — Conclusion 

We have investigated a simple one dimensional model for the swelling or the collapse of a closed 
polymer chain of size N, representing the dynamical evolution of a polymer in a 0-solvent that is 
rapidly changed into a good solvent (swelling) or a bad solvent (collapse). 

• In the case of swelling, the initial phase is to a large degree independent of the chain constraint. 
The density profile is parabolic and expands in space as t 1 ^ 3 . This regime is well described by a 
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Figure 12: Evidence for aging in the structure factor, for N = 10000, q = 5, T = 5, averaged over 10 
initial conditions. The curves for t w = 10 4 , 2 10 4 and 3 10 4 superimpose when plotted as a function 
of t/t^ 4 . The curve with t w = 1000, however, has a different shape. The mechanism proposed in the 
text suggests that the scaling in t/t!^ should indeed only be approximate. 



Flory-like continuum equation for the monomer density. The dynamics slows down after a time 
oc iV 2 when the chain becomes stretched, and the polymer gets stuck in metastable 'zig-zag' 
configurations (that we have enumerated), from which it escapes through thermal activation. 
The final stage of the stretching is therefore found to be logarithmic in time: the size of the 
polymer grows as Vint. 

• In the case of collapse, the chain very quickly (after a time of order unity) breaks up into small 
'pearls' or clusters of monomers. The evolution of the chain towards the fully collapsed states 
proceeds through a slow growth of the size of these metastable clusters, again evolving as the 
logarithm of time. We have enumerated the total number of metastable states as a function 
of the extension of the chain, and deduced from this computation that the radius of the chain 
should decrease as a double logarithm of time, i.e. l/ln(lnt). We have also found the total 
number of metastable states with a given value of the energy. We also obtained the distribution 
of cluster sizes, that we compared to simple one dimensional coalescence models. Finally, the 
aging properties of the dynamical structure factor was investigated numerically. The subaging 
behaviour that we found was attributed to the tail of the distribution at small cluster sizes, 
corresponding to anomalously fast clusters. We have argued that this mechanism for subaging 
might hold in other slowly coarsening systems. 

The system studied here falls in the category of non disordered, finite dimensional models exhibit- 
ing glassy dynamics. The list of such models keeps increasing with time p8|| , although the present 
model is quite realistic (no ad- hoc dynamical rules or long-ranged interaction). Actually the chain 
constraint does not play a very important role for the dynamics of collapse which is, as mentioned 
above, quite close to the one- dimensional version of the Backgammon model. 

In real polymeric systems, the presence of a chain would change quantitative features such as 
diffusion coefficients. Polymer melts have been shown to behave very similarly to structural glasses 
such as binary mixtures of Lennard- Jones particles |TIJ. These polymeric systems are out of equi- 



librium at long times and attempts have been made to describe them within the framework of the 
Mode Coupling Theory [ID, Yl\. 



A limitation of the present model in the collapsing regime is the unrealistic nature of the in- 
teraction potential: there is no hard-core constraint, and the range of the attractive part is finite. 
For realistic potentials, the total energy of a cluster containing n monomers will eventually grow 
like n for large n and not like n 2 . Correspondingly, the energy barrier to remove one monomer 
from a large cluster is finite, instead of being of order n. Furthermore, the long-ranged part of the 
attractive potential can speed up the dynamics by making the large clusters move towards each 
other, an effect absent in the present model. Nevertheless, the model used here could be justified for 
colloidal particles, of the type studied in fl3f . In this system, clusters of particles form and slowly 



coalesce, this leads to a very clear aging of the dynamical structure factor. Hints about the presence 
of 'micro-earthquakes', possibly due to the disappearance of large clusters, have been reported. A 
closer comparison between our model and such experiments would be very fruitful, but might require 



the inclusion of elastic deformations, as suggested in [13 
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Appendix A: A variational method for the swelling of a poly- 
mer chain 



In |§ a polymeric chain with given two-body interactions is approximated by a 'variational' gaussian 
chain with a time-dependent Kuhn length. In other words one starts from the original chain with 
hamiltonian H and its associated Langevin equation: 
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and then one tries to replace it with another system characterized by H v : 
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rj(s, t) is a gaussian noise with 



( V (s,t)r](s,t')) = 2D8(t-t') 



(49) 
(50) 

(51) 
(52) 

(53) 



and Tn = where D is the diffusion coefficient. 

u k B T 

In order to find the effective Kuhn length that would represent best the initial system, the radius 
of gyration is chosen to coincide in both models, to first order in x( s >t) = r ( s >^) — r^ v '(s,t) and 
W = H-H V . 

This leads to an implicit equation for a(t): 
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<r^(s,t) X (s,t) >=0 



(54) 
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We have used periodic boundary conditions, introducing uo n — y . 
The whole calculation is discussed in detail in |J, and we only reproduce the final equations 
which appear in the form 1 + 11 = 0, with, in the case of two-body interactions: 
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(55) 
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We have set: 



a(r, U) =4^(1- cos 2p7r*7) 
p>i 



a l -2Dj*f(r) , 

Nul N Jo 



dXe -2DJ>{f(T)-fW) 



(57) 



and /(f) = jj 



We are interested here in an intermediate time regime, and choose to keep the time dependence 
in a(t) such that R(t) = a(t)Nv <C and look for a power-law scaling: a(t) = N^t a . Then the 

reduced variable u = with 7 = ^ [ 2 (d+2) ~~ fA * s ^ ne re l evan t parameter in this time regime; we are 
interested in the situation where one can take the limit N —>■ oo, keeping u finite and much smaller 
than unity. 

A succession of approximations gives the following behaviours for the quantities / and II: 
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where the parameters used above are defined as: 
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2a, 



P = 7 + fi ~ 2 X 7, 
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(58) 



(59) 



Finally this allows to find the result announced in the main text for the case of repulsive interactions: 

a=i£j (d<2), (60) 

so that, in d = 1, -R(t) ~ t^A^Ti. 



Appendix B: Calculation of the number of metastable states 
in the ID Kawasaki Ising model 

Cornell, Kaski and Stinchcombe |y| have studied the dynamics of the one dimensional Ising model, 
using the Kawasaki prescription: the total magnetization is fixed and the only moves allowed are 
exchanges of two antiparallel nearest-neighbour spins. It is best suited to model binary alloys. At low 
temperature, the model has been shown to exhibit a glassy-type behaviour. At finite temperature, 
the system evolves in time by condensation of domains of parallel spins (if J is the exchange coupling 
constant, there is an energy gain of 4 J for a spin to join a domain by switching sites with its nearest 
neighbour); conversely, the evaporation of a spin from a domain costs 4 J. The diffusion of a spin 
performing a random walk in a domain consisting of spins of opposite sign costs no energy. 

At zero temperature, all evaporation processes are frozen and only condensation and diffusion of 
spins can occur. However, in order for two domains of same sign to coalesce, the evaporation of the 
boundary spins is necessary; as a consequence, at T — 0, many states are frozen. 

The frozen states are in fact very easy to characterize: they consist of alternating positive and 
negative domains, each of them containing at least two spins. Let M be the total number of positive 
domains; for i = 1, . . . , M, nf is the number of + spins in the positive domain i. Between two 
of these domains there exists a domain of negative magnetization. We label also these domains by 
i — 1, . . . , M, and n~ is the number of — spins in the negative domain i. 



Since iV is the total number of spins and Nb the total magnetization (the letter m being used in 
the text for the chain constraint), we have the two conservation relations: 

M M 

N = J2(nt + nry, Nb = J>+ - n~). (61) 

i=l i=l 

The number of frozen states for M fixed and fixed magnetization is given by the following expression: 

N N N N / M \ / M \ 

N F (M,b) = (Y, ... E)(E ••• E)MEk + + 0-JvU EK + -nn-JV6 ■ (62) 

n+=2 "+=2 71^=2 ™m= 2 

The calculation can be handled easily by exponentiating the (^-functions, and taking the N — > oo 
limit. 
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Then a saddle point approximation on = «A and z-i = ifi leads to a simple expression for the 

>M_ 

N 



entropy of frozen states at fixed magnetization per spin m and fixed fraction of domains 5 = 2 M- 



-^-\nN F (M,b) = S{b,5) = -{-l-b + 25)ln{l + b-25) + -{l+b-5)ln(l + b-5) 
./* 2 2 

+ -(-1 + b + 25) ln(l -b-25) + -(l-b-5) ln(l - 6 - 5) - 5\n5. 
2 2 

We then determine the optimum 5* for 5 at fixed magnetization b, in order to find the total 
entropy for a given b. This leads to an implicit equation: 

(l + b-25*) 2 (l-b-25*) 2 = 

5* 2 (l + b-S*)(l-b-S*) 1 } 

Case 6 = 

The complexity of the Kawasaki model is then easy to compute, 

S K (b = 0) = (-1 + 25) ln(l - 25) + (1 - 5) ln(l - 5) - 5 In 5, (64) 

The sum over 5 can be performed using a saddle-point approximation, leading to a total number 
of frozen states: 

(~<N 

N*(b = 0)~-=, (65) 



where the complexity per spin is C ~ 1.618. This is slightly less than the number of closed random 
walks for which C = 2: each metastable configuration in the Kawasaki problem can be seen as a 
random walk that goes at least twice in the same direction before changing direction. 



The case of zero magnetization is of special interest because it allows to give an estimate of the 
complexity of the polymeric chain in its stretching regime. As was pointed out in section 3.3, as the 
chain swells, it is trapped in metastable configurations, consisting of fully stretched segments of the 
chain going alternatively to the left or to the right. It is easy to draw a correspondence between a 
segment going to the right (respectively, to the left) and a domain of plus (respectively, minus) spins. 
The constraint that the chain is closed requires that the number of steps to the right is equal to the 
number of steps to the left, which is equivalent, in the spin language, to b = 0. 

In order to count the total number of metastable states for the chain, one should actually require 
that each stretched segment must contain at least one monomer, so that the sums in ( |6"2]) should 
start at nf = 1. However, there must be an energetic bias in order for each end of segment to be 
stuck at its position, namely that the number of monomers on its 'inward' neighbouring site along 
the chain should be more important. As explained in the main text, we assume that this is true with 
probability 1/2, and this adds a multiplicative factor in formula (|62|) . Finally the complexity of 
the chain is: 

S = -<51n2 + (5- l)ln(l -5) -5\n5, (66) 

This relation is illustrated on Fig. 4 as a function of R, given that R ~ Ms. The sum over 
5 can be performed using a saddle-point approximation, leading to a total number of frozen states 
N F ~ C N /VN, where C = 3/2. 

To summarize, the calculation reveals an interesting analogy: the number of states of a chain 
containing at least k monomers on each stretched segment is the same as the number of the frozen 
states in the ID Kawasaki spin model at zero magnetization and zero temperature (with all domains 
containing at least k spins), and is in turn equal to the number of one- dimensional closed random 
walks biased in such a way that one has to walk at least k times in the same direction before changing 
direction. 

In order to know the number of metastable states of the chain with at least k monomers on each 
segment, one must add an energetic bias introduced via the additional factor ^ir in equation (p2]). 
To complete the analogy let us formulate two remarks. 

First, the number of closed random walks in one dimension (2 N / y/N) is recovered when one uses 
the above formula (|62| ) with sums starting at nf = 1. 

Second, one can also vary the degree of bias k of the random walk by starting the sums in ( p2|) at 
nf = k. Physically this corresponds to metastable states over a time scale t such that \nt = (3vk. If 
in addition, we incorporate the energetic bias factor ^it, this allows to follow with time the number of 
metastable states at finite temperature as the system coarsens, and relate it to the radius of gyration 
of the chain at long times. 

We obtain the following entropy of metastable states, as a function of the parameter u = Tint: 

S(u) = -5 In 2 + (-1 + u5) ln(l - u5) + (1 - (it - 1)5) ln(l -(u- 1)5) -5\n5, (67) 

This entropy is maximum at a certain value of 5, which, to leading order in u (for u large), is 
given by 5* ~ -. We expect that the number of zig-zags in the polymer will adjust at this value, 
corresponding to one of its most probable configurations. As was explained in the main text, R and 
5* are then related by 5 ~ N/R 2 , so that R ~ VN In t. 

Finally, let us mention that the whole calculation is in spirit (and also in the analytical methods 
involved) very similar to the calculation of the number of metastable states of a collapsing chain, as a 
function of the degree of freezing k* = T In £; section 4 of this article is devoted to such a description. 



Appendix C: Number of collapsed metastable states as a 
function of the energy 

We start by the expression of the number of A;*-stable configurations with M clusters and iV 
monomers, at fixed energy E, given by: 

N N / M \ / M i \ 

M p , h .{E,M\N)= ••• E S[J2(ni-2) + 2L-N)5[J2-i,n i (n i -l)-(L-M)-E). 

ni =k* n M =k* \i=l ) \i=l J 

(68) 

In the following, we set e = -j| < and limit ourselves to the case m—1. 

Let us first note that for a given k*, the maximum energy per particle is equal to e max = — + 
1) + y. It is reached for a configuration consisting of M clusters of k* monomers each, each of 
energy —k*(k* — l)/2 and the remaining L — M sites are occupied by 2 monomers, each of these 
sites contributing with an energy equal to —1. The conservation relation N = Mk* + 2(L — M) 
leads to 5 max = M ™ ax = (1 — 2/p)/{k* — 2), and the energy is calculated accordingly, leading to 
e max = —\{k* + 1) + — . In the following, we will only need to consider values of e in the interval 

] O^- 1 , C ma x] • 

From the expression of J\f p ^*(E, M\N) above, one may use the integral representation of the 
(^-functions to find that: 

Af Plk .(E, M\N) = dX J dYe- NF{x ' Y) . (69) 

with 

F(X,Y)=X(l--)+Y(e + -)-6hiE(X,Y), (70) 
P P 

and 

oo 

E(X,Y)= e x "-^«(™+3). (71) 

n=k*-2 

Our procedure is to find an estimate of this quantity via a saddle-point approximation, so that 
the entropy of metastable states at fixed 8 is given by: 

where the last step is to maximize this expression with respect to 5 so that the final result is 

S p (e) = S p (5*,e). 



Limit of low energies e — > — oo 

In this part, it is more convenient to perform a change in variables Z = X — 3Y/2 and Z' = —Y/2. 
Moreover, numerical investigation suggests that for |e| large, Z* should be small and Z'* of the form 
Z'* ~ — aZ* 2 , a > 0; so we set in the following Z' ~ — aZ 2 , a > 0, so Z and a are the new variables: 



F(Z, a) = Z(l --) + aZ 2 (3 + 2e--)-S In S(Z, a), 
P P 



(73) 



and: 

oo 

E(Z,a)= Yl e~ aZ2n2+Zn (74) 



n=k*-2 



We then approximate the discrete sum E(Z, a) by an integral, which is valid if we find a posteriori 
that Z is small: 

! poo 

E(Z,a) ~ --G(a), G(a) = due'™ ~ u (75) 
Z Jo 

Up to this order in 1/Z the results are independent of k*. The combined saddle-point equations 
dF/dZ = and dF /da = lead to the simple relation: 

*.(!-?) = (76) 

Continuing the algebra with the assumption that a* 7^ does lead to inconsistent values for 5* 
(5* > \). We conclude that the saddle-point solution, if there is one, is stuck at a* = 0. The rest of 
the calculation is then straightforward: we obtain Z* = —6/(1 — -), and since G(0) = 1, 



Finally, in order to determine 5*, we compute: 

^ 2 P (i - M 



^^ = ln(1 -p } - ln 



(i-25 P y 



(78) 



In the limit e — > —00, the interval allowed for 5 is [0; 5 max = —(1 — - p ) 2 /2e\. In this interval, S P (S, e) 
is an increasing function of 5, and therefore reaches its maximum at the upper bound. Therefore: 
5* = <5max- Note that this implies that Z* — (1 — |)/2e is indeed small as was assumed in this 
approach. We get the final result for the tail of the entropy of metastable states : 

S P (e) -^(1 - -) 2 + -(1 - -) 2 [-l + \ ln(£(l - -) 3 / 2 )] + o(±) (79) 



pep 2 4 p e 2 



Limit of high energies e — > e 



max- 



We start again from the expression for the free energy in terms of the variables X and Y. We assume 
that Y* to be large, and a good approximation is to keep only the first two terms in the sum S(X, Y), 
corresponding to n = k* — 2 and n — k* — 1. Then, 

F(X, y) = X(l - - - 5(k* - 2)) + y(e + - + -<?(Jfe* - 2)(F + 1)) - 51n(l + e x - fc * y ). (80) 
P P 2 

The saddle-point equations dF/dX = and dF/dY = lead to a simple system for x = e x and 

Y 

y = e r : 



x 1 - I - - 2) x €+ l + l 5 (k*-2)(k* + l) 



,k* 



5(k* - 1) - 1 + 2- ' e + i + i<J(ib* - l)(k* + 2) ' 



(81) 



From this system we obtain a relation between X and Y: 



5(k* - 1) - 1 
as well as a consistency relation between 5 and e: 

(k* - l)(k* - 2) 



1 - 2 - <y(ife* - 2) 

X = FY + ln| — -g - 2 ; 1 , (82) 



(83) 




Let us note that at the upper boundary of the interval studied, i.e. for e = e max — — ~(k* + 1) + — , 
this consistency relation gives S = S max = (1 — 2/p)/(k* — 2), which is indeed the exact result. Then 
the expression for entropy can be computed in a simple way: 

2 fi-l-S*(k*-2)\ ( 5* 

SJe) = -(1 <T(fc*-2))ln - ./ -— T}+S*\n\— - 

pw v p v 77 I S*(k* - 1) - 1 + - I \6*(k*-l)-l 

" g ln { (i-25yy ) + -p ln {T^W-p) (84) 

This gives an estimate for the entropy of metastable states near e max (k*) for a given k*. Note that 
the expression found above is actually defined only in the interval [e max {k* + 1) = —\{k* + 2) + 
^p-; £max(k*) = —\{k* + 1) + — ], beyond which the present approximation breaks down. 

To summarize, for each k*, or equivalently for each u = Tint, we found the behaviour of the 
entropy of metastable states both for very low energies (e — > — oo), where it is independent of k*, and 
in the vicinity of the maximum energy e max (k*). Thus the full curve for S p (e, k*) is an extrapolation 
between these two asymptotic regimes, as shown in Fig. 10. 
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